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Abstract 

O 

The relativistic Hartree-BCS theory is applied to study the temperature 

o ' 

dependence of nuclear shape and pairing gap for 166 Er and 170 Er. For both 
the nuclei, we find that as temperature increases the pairing gap vanishes 
^ ' leading to phase transition from superfluid to normal phase as is observed in 

c ' 

nonrelativistic calculation. The deformation evolves from prolate shapes to 
• <— i . 

spherical shapes at T ~ 2.7 MeV. Comparison of our results for heat capacity 



with the ones obtained in the non-relativistic mean field framework indicates 
that in the relativistic mean field theory the shape transition occurs at a 
temperature about 0.9 MeV higher and is relatively weaker. The effect of 
thermal shape fluctuations on the temperature dependence of deformation 
is also studied. Relevant results for the level density parameter are further 
presented. 

PACS numbers: 21.10.Ma, 21.60.-n, 27.70.+q 
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I. INTRODUCTION 



The relativistic mean field (RMF) theory (l|-|3j has been very successful in describing 
the ground state (zero temperature) properties of nuclei over the entire periodic table in 
recent years. The binding energies, charge radii and the ground state deformations are 
reproduced very well; the charge distributions also compare extremely well with the exper- 
imental data. This theory has proved to be very fruitful in explaining various details 
of exotic nuclei near the drip lines. In contrast to the nonrelativistic models, the RMF 
theory uses a single set of parameters to explain all these properties. To our knowledge, 
this approach has not yet been exploited to understand the properties of hot nuclei except 
some preliminary investigations for closed-shell nuclei 0. The response of nuclear shapes 
to thermal excitations, for example, has been experimentally studied from the shapes of 
the giant dipole resonances (GDR) built on excited states ||[|. Theoretically such finite 
temperature effects have been studied till date in the nonrelativistic framework such as the 
Hartree-Fock-Bogoliubov (HFB) theory flT0| , |TT| and the Landau theory of phase transition 



12, pjfl. Such theories qualitatively explain, for example, the temperature evolution of nu- 



clear shapes. Recent experiments indicate, however, that a quantitative estimate of the 
persistence of the ground state deformation ]14[ with temperature may be missing in some 
cases. 

Against this backdrop, we undertake the study of the thermal response to nuclear prop- 
erties, particularly the deformation and the level density in the RMF framework. The 
HFB calculations employ a model hamiltonian in a limited model space with pairing-plus- 
quadrupole interaction. This may not be very realistic once temperature comes into play. 
Moreover, to make the calculations numerically tractable, an inert core is assumed. This 
may be questionable at moderately high temperatures. The RMF approach we use is effec- 
tively free from these limitations. The model space employed is sufficiently large and all the 
nucleons are treated on equal footing. 

In the present work, we employ the nonlinear a — uj — p version of the RMF theory [H . 
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In absence of a simple relativistic prescription for the pairing interaction, it is introduced 
phenomenologically. We take two representative systems, namely, im Er and 170 Er. We in- 
vestigate the thermal evolution of the nuclear shapes and the pairing gaps. The temperature 
dependence of the specific heat as a possible signature of phase transition in pairing and 
nuclear shapes is explored. The temperature variation of the nuclear level density parameter 
which has a pre-eminent role in understanding nuclear reactions is further studied. Effects 
of thermal fluctuations of the nuclear shapes on the deformation and the nature of the phase 
transition is also discussed. 

The theoretical framework used is briefly discussed in Sec. II. The results and discussions 
are presented in Sec. Ill and the Sec. IV contains the concluding remarks. 



II. FORMALISM 



The Lagrangian density for the nucleon-meson many-body system |2|] is taken as 
C=% (i-fdn -M)% + ^ad^a - U(a) - g a %(r% 

-gfii-fp^i - \f^F^ - e^S-^^j. (1) 

The meson fields included are the isoscalar a meson, the isoscalar-vector uj meson and 
the isovector-vector p meson. The arrows in eq. ([TJ) denote the isovector quantities. The 
Largrangian contains a nonlinear scalar self-interaction term U(a) of the a meson: 

U(a) = l -m 2 a a 2 + ^g 2 a 3 + (2) 



This term is important for appropriate description of the surface properties | 15|| . The quan- 
tities M, m a , m^j and m p are the nucleon, a, to and the p— meson masses, respectively, while 
9a, 9ui, 9p and e 2 /4iT = 1/137 are the corresponding coupling constants for the mesons and 
the photon. The field tensors of the vector mesons and of the electromagnetic fields have 
the following structure: 



W = d^uf - d u u», (3) 
RT = - d v (T - g p (jT x (T), (4) 
= d»A v - d v A». (5) 

The variational principle gives the equations of motion. The mean field approximation is 
introduced at this stage by treating the fields as c— numbers or classical fields. This results in 
a set of coupled equations, namely the Dirac equation with potential terms for the nucleons 
and the Klein-Gordon type equations with sources for the mesons and the photon. For the 
static case, along with the time reversal invariance and charge conservation the equations 
get simplified. The resulting equations, known as RMF equations, have the following form. 
The Dirac equation for the nucleon is 

{-ia • V + V(t) +{3[M + S(r)]}% = e^, (6) 

where V(r) represents the vector potential 

V(r) = ^o(r) + g p T 3 p (v) + e iI_^U (r), (7) 

and S(r) is the scalar potential 

S(r)= 9a a(r), (8) 

which contributes to the effective mass as 

M*(r) =M + S(r). (9) 

The Klein-Gordon equations for the mesons and the electromagnetic fields with the 
nucleon densities as sources are 

{-A + ml] a(r) = -g aPs {r) - g 2 a 2 {v) - g 3 a 3 (r), (10) 

{-A + ml}uo(T)=g u p v (r), (11) 

{-A + m 2 p }p (r)=g pP3 (r), (12) 

- AA (r) =ep c (r). (13) 



The corresponding densities are 

p s = ^mViVi, 

i 
i 
i 

Pc = T,**\^-^-*i- ( 14 ) 

Here the sums are taken over the particle states only, i.e., the negative-energy states are 
neglected. The partial occupancies (n,) at finite temperature in the constant pairing gap 
approximation (BCS) is 

1 
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1-^(1-2/(6^)) 



(15) 



with f(£i,T) = 1/(1 + e €4//T ); q = \J — A) 2 + A 2 is the qusiparticle energy where e« is 
the single-particle energy for the state i. The chemical potential A for protons (neutrons) is 
obtained from the requirement 

*Em = Z(N) (16) 

i 

The sum is taken over proton (neutron) states. The gap parameter A is obtained by min- 
imising the free energy 

F = E-TS, (17) 

where 

E(T) = J2 ^ + E « + e °nl + E W + E P + E C + E pair + E CM - AM, (18) 

i 

and 

S = - E {/Mi + (1 - /i)!n(l - fi)] , (19) 

i 

with 



E a = -~g a Jd 3 rp s (r)a(r), (20) 
Knl = ~\J d\^-g 2 a\v) + \g z a\r)} , (21) 
E u = -\g u [ d 3 rp v (r)uj (r) } (22) 



2' 



E P 



E, 



c 



-^p/^P3(r)p°(r), (23) 
-f ! d 3 rp c (r)A°(r), (24) 

<37T J 



Epair — q i (25) 

E CM = -hco = -hlA-V\ (26) 
Here G and A are the pairing strength and the mass number respectively. The single-particle 



energies and the fields appearing in eqs. (|T8|) - ( p3|) are obtained from the self-consistent 



solution of eqs. @ - (1T2f) . We generate these self-consistent solutions using the well tested 
basis expansion method as described in detail in Refs. P,|16 . 



The self-consistent solutions for the RMF equations at a given temperature is obtained 

by minimizing the free energy F which yields the equilibrium or the most probable value 

of the quadrupole deformation (/J^) and the proton (neutron) pairing gaps A p (A n ). To 

incorporate the effects of the thermal fluctuations, a constrained calculation is performed to 

generate the free energies for the deformations away from the equilibrium value. Average 

value of a quantity then can be obtained as 

J d(3 2 Q(f3 2 )e-^y T 
f d/3 2 e- AF (fc)/T 

where AF = F{(3 2 ) — F(f3 2 ) and 0(f3 2 ) is the expectation value of the operator O at a fixed 
value of (3 2 and T. The quantity e _AF//T is a measure of the probablity for the nucleus to 
have deformation f3 2 at the temperature T. 



III. RESULTS AND DISCUSSIONS 



We have chosen 166 Er and 170 Er as two representative systems for our calculations. 
The results presented are obtained using the NLSH parameter set [|16| for the values of the 
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coupling constants and masses for the mesons and nucleons. The pairing strength G is taken 
to be 29/A and 21/ A for protons and neutrons, respectively; they reasonably reproduce the 
observed pairing gaps at zero temperature obtained from odd-even mass differences. The 
values of the chemical potential and the pairing gap are determined using all the single 
particle states up to 2%Uq above the Fermi surface without assuming any core, i.e., the mean 
field solution is generated by taking into account all the nucleons in the systems considered. 

At finite temperatures, the nucleus is not strictly in thermal equilibrium. In order to 
treat the system as an isolated one in equilibrium, one has to take into account corrections 
due to excitations of nucleons in the continuum. This is generally done through a subtrac- 
tion procedure, by treating the liquid plus vapour phase together and then the vapour phase 
seperately |17| |. However, calculations of the nuclear level density parameter at finite tem- 



perature show that the results are insensitive |18j to the continuum corrections for T upto 
~ 3 MeV; we have therefore not included the effects due to the continuum in the present 
calculations as temperatures above it are not relevant. 

In subsec.3.1, the results for the most probable values (i.e, the mean field values ) of the 
different observables are presented. With increasing temperature, thermal Auctions build 
up which may shift the average values from the most probable ones. In subsec.3.2, we 
discuss the results with the inclusion of thermal fluctuations. Fluctuations in both the fa 
and 7- degrees of freedom of the nuclear shape should be incorporated, however, because 
of simplicity and computational economy, we have included only the /3 2 - fluctuations in the 
present calculations. 



A. Mean field results 

We have calculated the most probable values of the quadrupole deformation parameter, 
neutron and proton pairing gaps, heat capacity and the level density parameters for the 
aforesaid systems function of temperature. 

In Fig. 1, the variation of the quadrupole deformation parameter (5 2 as a function of 
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temperature is shown. The ground state deformations obtained here are almost close to 
the ones calculated in Refs. [|K],|n|] using nonrelativistic mean field theory with pairing plus 
quadrupole (P+Q) interaction. The nature of the temperature dependence of the deforma- 
tion in the two cases are also not very different. However, whereas in the nonrelativistic case 
the phase transition from prolate to spherical shape occurs at T ~1.8 MeV, in the present 
case the said transition is found at a higher temperature, T ~2.7 MeV. The signature of the 
phase-transition can be inferred by examining other observables like heat capacity which we 
consider later. 

In Fig. 2, the results for the evolution of the neutron and proton pairing gaps with 
temperature are displayed. For both the nuclei the pairing gaps monotonically decrease 
with the increase in temperature. It is seen that for both the nuclei, the neutron pairing 
gaps vanish almost at the same temperature T ~0.4 MeV; similarly the proton pairing gaps 
vanish at T ~0.45 MeV. These results are not much in variance with those obtained in the 
nonrelativistic P+Q model. 

The vanishing of the nuclear deformation and the pairing gaps with temperature indicates 
that there is a shape transition from prolate to spherical and also a transition from the 
superfluid phase to the normal phase. To understand the nature of the transitions, the heat 
capacity at various temperatures are calculated. The specific heat C for a given temperature 
is obtained using 

C{T) = °— (28) 

where E* is the excitation energy of the nucleus. In Fig. 3, the temperature variation of the 
specific heat for im Er and l70 Er nuclei is plotted. At T ~ 0.4 MeV, two closely seperated 
peaks in the C(T) curves are seen. They correspond to the dissolution of the neutron and 
proton pairing gaps. These are the characteristic signatures of second order phase transitions 
from superfluid to normal phase. It is further noted that the heat capacities for the two 
systems have broad bumps at T ~ 2.7 MeV. This signifies a weak second order phase 
transition corresponding to the transition of the nuclear shape. This result is at variance 
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with those obtained in the P+Q model or in the nonrelativistic calculation |T2[ based on the 
Landau theory of phase transitions where a strong second order phase transition is observed. 

We now present the results for the temperature dependent level density parameter a. 
The parameter a can be obtained using the excitation energy and the entropy as follows, 

E* = aT 2 , (29) 
S = 2aT. (30) 

The parameter a obtained using above eqs. (|29"D and ([30|) would be equal provided it is 
independent of temperature |fL9"| . In Fig. 4, the inverse level density parameter K = A/a (A 
is the mass number of the nucleus) is plotted as a function of temperature. The subscripts 
E and S are used to distinguish the two definitions given by eqs. ( p9[) and (|30"D , respectively. 
At lower temperature, both Ke and Ks shoot up due to the collapse of the pairing gaps. At 
higher temperatures, Ke and Ks are quite close to each other and there is no appreciable 
variations in their values in the temperature range considered. This is due to the weak 
transition of the nuclear shapes from deformed to spherical at T ~ 2.7 MeV. 

In order to test the sensitivity of the results presented above to the choice of the parameter 
set and the model space, we have repeated the calculations for the nucleus im Er with the 



NL3 parameter set pO|,^TJ. In comparison to the results for NLSH parameter set it is 
found that the NL3 parameter set gives slightly larger (~ 10%) value of the ground state 
deformation. For T > 1.5 MeV, the values of deformation obtained in both NLSH and NL3 
parameter sets are almost identical. The temperature evolution of the pairing gaps and level 
density parameters are practically the same with the two parameter sets. We extend the 
model space to include single particle states upto 3huo above the Fermi surface instead of 
2hu Q as used above. For this extended model space the pairing strength G is adjusted to 
reproduce the ground state pairing gap. We do not find significant changes in the values of 
deformation parameter and the pairing energy in the temperature range of interest (T ~ 3 
MeV). In the extended model space, the value of the inverse level density parameter is also 
nearly the same. The model space used is thus found to be sufficient for our calculations. To 



estimate the importance of the continuum corrections to the inverse level density parameter 
we calculated the occupancy, n^ + \ of the single-particle states with positive energy. For 
T < 1 MeV practically there is no particle in the positive energy states (n^ = 0) and at 
the highest temperature of interest studied n^/A = 0.018 which is very small. Thus it is 
expected that the continuum corrections may not play an important role in the temperature 
range we study. The continuum effects may grow stronger for T > 3 MeV, however, this is 
beyond the transition temperatures and so we have not taken this into account. 

The shape transition temperature and the inverse level density parameter obtained in 



the present model are higher compared to those calculated in the (P+Q) model |[T0| , [TT 
Moreover, in the later model, the deformation falls to zero sharply whereas in the present 
case, this fall is comparatively a little slower. The origin of these differences can be attributed 
to the single-particle level spectra. In Fig. 5, the energy levels for the different proton orbitals 
around the Fermi surface for the nucleus 170 Er at T=3.0 MeV are displayed. Also the proton 
single-particle energies in the (P+Q) model are shown at this temperature. In both models 
the orbits are spherical at this temperature. The striking difference between the orbitals in 
the two models is the large energy gap across the Fermi surface in the RMF model. Similar 
is the case for neutrons. These tend to reduce the entropy and the excitation energy and 
hence the level density parameter. The increase in the shape transition temperature (by ~ 
0.9 MeV) in the RMF model as compared to the P+Q model may be traced back to the 
specific level structure and the gap across the Fermi surface in the two models and their 
thermal response. 

B. Results with thermal fluctuations 



Experimental data on giant dipole resonances built on excited states give estimates of 
the average values of nuclear deformation at finite temperatures. They can be calculated 
by taking into account the thermal fluctuation effects around the 'most probable' value. 
Nuclei being finite systems, one expects that thermal fluctuations would play an important 
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role in the quantitative estimation of various nuclear observables at finite temperatures. In 
fact, it has been shown in Ref |22[ that the experimental data on GDR built on the excited 
states can be explained reasonably only if the effects from the thermal fluctuations of all the 
quadrupole degrees of freedom for nuclear shape are included. In the following, we consider 
the thermal fluctuations only in the /?2 degree of freedom. Inclusion of fluctuations in the 
7-degrees of freedom is very involved and computer-intensive and therefore not considered 
presently. We perform constrained calculations for free energy at fixed deformations 02- 
The free energy surface shows multiple minima for temperatures below the shape transition 
temperature. The average value of fa is calculated using eq.(27). In Fig. 6, we have displayed 
the plots of free energy surface AF for l70 Er at a few temperatures ranging from 2.0 to 3.0 
MeV where AF is measured from the lowest minimum (prolate in the present case) at 
the corresponding temperature. For T = 2.0 MeV, the free energy surface has minima at 
/?2 = —0.16 and 0.23. The free energy difference between these two minima is 0.70 MeV (i.e. 
AF/T = 0.35). As the temperature approaches the shape transition temperature (T ~2.7 
MeV), both the minima in the free energy surface merge and flatten the bottom part of the 
free energy surface. This is evident from the middle panel of Fig.6 corresponding to T = 2.65 
MeV. With further increase in temperature, the free energy surface with single minimum 
at zero deformation broadens. The variation of the average deformation (^2) as a function 
of temperature obtained through eq . (|27|) for 170 Er is shown in Fig. 7. The average value of 
@2 differs significantly from its most probable value (see Fig. 1) for temperatures above 1 
MeV. The reason for such differences can be understood as follows. The free energy surface 
for 170 Er at low temperatures (T < 2 MeV) shows two well seperated minima. One of these 
minima lies at prolate deformation (F p ) and the other one is at oblate deformation (F a ). 
Their difference AF° = F Q — F p as a function of temperature is shown in Fig. 8. For T > 0.5 
MeV, AF° decreases monotonically to zero. The initial rise of AF° can be attributed to 
the quenching of the pairing correlations; it is seen that without pairing, AF° decreases 
uniformly with temperature. The decrease in AF° enhances the probability of finding the 
nucleus with oblate deformation and hence the average deformation is smaller than the 
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most probable one for temperatures below the transition temperature. The deviation of 
the average deformation from the most probable value is also governed by A(3 2 — P2 ~ P%- 
Here, (3 2 and (3 2 are the deformations corresponding to the prolate and the oblate minima. 
With rise in temperature, they come closer and ultimately merge as shown in Fig. 9. The 
persistence of j3 2 beyond the transition temperature is attributed to the asymmetry in the 
free energy surface as seen in the bottom panel of Fig. 6. A similar behavior for im Er is also 
seen. 

In the mean field approximation, the thermal evolution of the most probable value of the 
deformation is seen to be associated with two phase transitions, the first being a superfluid 
phase to normal phase transition at T ~ 0.4 MeV and the other a weak second order 
phase transition in the shape at T ~ 2.7 MeV. It would be interesting to see how thermal 
fluctuations affect these phase transitions. In Fig. 10, the specific heat for the system 170 Er is 
displayed as a function of temperature with inclusion of fluctuations. At low temperatures for 
T < 1.0 MeV, there is no perceptible change in the specific heat and thus the transition from 
superfluid to normal phase is not affected by fluctuations. At a little higher temperature, 
however, the broad bump seen in Fig.3 at ~ 2.7 MeV is seen to be shifted to ~ 1.6 MeV and 
becomes also a little sharper. The phase transition in nuclear shape thus occurs at a lower 
temperature and is not as weak as seen in the mean field approximation. Normally, one 
would expect the phase transition to get diluted due to presence of fluctuations. Here one 
sees the opposite. The evolution of the free energy and the energy profiles with temperature 
showing two minima (one prolate and the other oblate) in deformation space is responsible 
for this observation. Here the role of thermal fluctuations becomes significant only in the 
intermediate temperature domain (temperature in the range 1 - 2.5 MeV, approximately) due 
to decrease in AF° and Af3 2 (see Figs. 8 and 9). The delicate interplay of these quantities 
in the averaging augurs a change in the heat capacity signalling a somewhat sharper shape 
transition at a lower temperature compared to the most probable one. 

The effect of thermal fluctuations of the shape coordinate f3 2 on the inverse level density 
parameter K for the system 170 Er is shown in Fig. 11. The full line and the dashed line 
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correspond to results of our calculations for Kg and Ke-, respectively. The corresponding 
lines with circles represent those with inclusion of the fluctuations. It is seen that the 
fluctuations shift the minima for Kg and Ke from 2.7 MeV to 1.6 MeV, the corresponding 
transition temperature. Beyond the transition temperature, the upward slope for K is also 
found to be higher with fluctuations included. Similar is the behavior for im Er. 

IV. CONCLUSIONS 

The relativistic mean-field approach together with pairing effects included through the 
constant gap approximation has been applied to study some properties of axially deformed 
im Er and 170 Er nuclei at finite temperature. The temperature dependence of the quadrupole 
shape for both the nuclei are found to be practically the same; similar is the case for the 
pairing gaps. The deformation remains close to its ground state value (f3 2 ~ 0.3) for T less 
than 1.5 MeV and falls to zero quite sharply at T ~2.7 MeV. Pairing gaps vanish at T =0.4 
- 0.5 MeV leading to transition from a superfluid phase to normal phase. Compared with 
the results obtained in the nonrelativistic mean field framework, the relativistic calculations 
yield vanishing pairing gaps at practically the same temperature and the fall is a bit slower. 
The evolution of the heat capacity with temperature and the transition from superfluid to 
normal phase do not behave very differently. The nature of the shape transition obtained in 
the two approaches are also not too different except that in the present case, the approach 
of the deformation to zero is relatively slower. This slower fall results in a broad bump in 
the specific heat around the shape transition temperature (T ~ 2.7 MeV). 

Fluctuations are expected to influence the phase transitions. To explore these aspects, 
we have taken into consideration thermal fluctuations in the (3 2 — degrees of freedom. The 
effects of fluctuations are found to be imperceptible below T ~ 1 MeV and therefore the 
pairing transitions are not affected. However, the influence on the shape transitions is found 
to be very significant; the transition temperature drops down from ~ 2.7 MeV to ~ 1.6 MeV 
with inclusion of thermal fluctuations. It also becomes a little sharper. The specific nature 
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of the free energy surface in the deformation space and its evolution with temperature are 
responsible for such a behavior. Fluctuations in the 7— degrees of freedom would also have 
a role to play in the phase transitions mentioned. We have not included these presently as 
they are computationally too intensive. 

The authors gratefully acknowledge Prof. P. Ring for providing them with the computer 
code to generate RMF solutions for axially deformed nuclei at zero-temperature. 
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Figure Captions: 

Fig. 1: Variation of the most probable quadrupole deformation j3 2 as a function of temperature 
for 166 Er and 170 Er. 

Fig. 2: Temperature evolution of neutron and proton pairing gaps for im Er and 170 Er. 

Fig. 3: Variation of specific heat as a function of temperature for im Er and 170 Er. 

Fig. 4: Temperature dependence of inverse level density parameters K E and K s determined 
from the excitation energy and the entropy, respectively. 

Fig. 5: Spectra of the proton single-particle energies around the Fermi surface for 170 Er at 
T = 3.0 MeV in the RMF and the (P+Q) models. The single-particle energies e are 
shown relative to the respective chemical potentials. 

Fig. 6: Plots for the free energy surface at various temperatures for 170 £V. The free energies 
at different temperatures are measured from the respective lowest minimum. 

Fig. 7: Average value of quadrupole deformation as a function of temperature for 170 Er. 

Fig. 8: The free energy difference between the prolate and the oblate minima as a function of 
temperature for 170 Er. 

Fig. 9: The difference A[3 2 between the deformations at the prolate and the oblate minima as 
a function of temperature for 170 Er. 

Fig. 10: Variation of average specific heat as a function of temperature for 170 Er. 

Fig. 11: The inverse level density parameter for the system 170 Er. The full and the dashed 
lines correspond to Kg and Ke, respectively. The lines with filled circles represent the 
corresponding results with inclusion of thermal fluctuations. 
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